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■ Propagators of the diagonal and the off-diagonal gluons are studied 

■ numerically in the Maximal Abelian gauge of SU{2) lattice gauge theory. 
^ . It is found that in the infrared region the propagator of the diagonal 
D ! gluon is strongly enhanced in comparison with the off-diagonal one. The 

', enhancement factor is about 50 at our smallest momentum 325 MeV. We 

' have also applied various fits to the propagator formfactors. 

5^ ; 1 Introduction 

The propagators of the fundamental fields play important role in the understanding of 
the physical structure of any quantum field theory. The gluon propagator in QCD is well 
known in perturbation theory, i.e. at large momenta. On the other hand its form in the 
infrared region has not been fixed so far although it has been intensively studied both 
analytically and numerically using lattice regularization. Analytical results range from 
the infrared divergent IJ to infrared vanishing [2] propagator (see also recent review 
The recent lattice investigations jH El 13 E] excluded the infrared divergent behavior 
leaving open the possibility of the infrared vanishing propagator. Another recent study 
[Oj - made in the Laplacian gauge which is free of Gribov copies - provided some support 
to the form of the propagator with dynamically generated gauge invariant mass proposed 
in dni. 

It is widely believed that the knowledge of the infrared behavior of the gluon propa- 
gator is crucial for understanding of the confinement problem. At present there are two 
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competing scenarios of confinement: condensation of monopoles ^I] or center vortices jT2]- 
The Maximally Abelian (MA) gauge is the most convenient for demonstration of the dual 
superconductor nature of the gluodynamics vacuum (see, e.g., jTS| for a review). The first 
study of the MA gauge gluon propagator in the coordinate space was made in [T3|. It 
was found that the propagator of the off-diagonal gluons is exponentially suppressed at 
large distances by the effective mass about 1.2 GeV. Thus the findings of 14 support the 
Abelian dominance in gluodynamics P^ll6j . The mass gap generation for the off-diagonal 
gluons was further studied analytically in Ref . [TTl HH] . 

In this paper we consider the propagators in the momentum space which allows a 
detailed investigation of their infrared properties compared to the coordinate space. Our 
preliminary results were published in [19]. We describe the gauge fixing and present defi- 
nition of propagators in Sections El and El Section (His devoted to discussion of numerical 
results and the last Section contains our conclusions. In Appendix we discuss the quality 
of the gauge fixing procedure, Gribov copies and finite volume effects. 

2 Gauge Fixing 

We use the standard parameterization of SU{2) link matrices Uu = cosife^^ and U12 = 
sin if e^^. The gauge fields are defined as follows 

where a"" are the Pauli matrices. In terms of the link angles one gets^: 

^^li^) = smyD^{x)cosXf,{x) , (1) 
X 

-Aj(x) = cosv9^(x) sin6'^(x) . 

We call v4^(x) the diagonal gluon field, and A^(a;), i = 1,2, the off-diagonal gluon field. 
The MA gauge condition in a differential form is pUj : 

^^T^Al{x) 

Note that here and below we are using the same notations for lattice and continuum fields. 
A nonperturbative fixing of this gauge amounts to the minimization of the functional 

FuAoiA] = J d'x{[Al{x)r + [Alix)]'} , 

which has the following lattice counterpart: 

^mag[^] = E^o^2^/^(^)- (3) 

X,fl 

^Note that in Ref. jHl the definition of the field A differs from Eq. (Q) by the factor of 2. 
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The MA gauge condition Q leaves U{1) degrees of freedom unfixed. To complete the 
gauge fixing we use a [/(I) Landau gauge. In continuum the Landau gauge condition is 

d,Alix) = 0. (4) 

In previous lattice studies [T3l[TH] to fix f/(l) Landau gauge the following lattice functional 

X,fJ, 

was maximized with respect to U{1) gauge transformations, 0^{x) —>■ 6^{x) + d^uj{x). 

In this paper we implement the following generalization of the U{1) gauge fixing func- 
tional © 

^and[^, V] = Y1 COS COS ^^(a;) , (6) 

X,fl 

which is consistent with the definition of in (Q). Contrary to the definition in Eq.(j3)) 
this condition implies that is transverse for any lattice spacing. In the continuum limit 
definitions (0) and © coincide with each other. 



3 Propagators 

We calculate the diagonal propagator 



Df:'ip) = {At{k)Ati-k)) , (7) 



and the off-diagonal propagator 

= {A^,ik)A-i-k)), (8) 
where the Fourier transformed field, A^^{k), is defined as follows: 

-^U^) = ;^E^"^'"^"^'^^;(^)' ^^='^ir' n, = 0,...,L,-l. (9) 



X 



The standard variables are 

2 . ak 



Pi, = - sm 



a 2 ' 

in terms of which lattice propagator of a free massive scalar particle in momentum space 
has a familiar form, D{p) oc + m^). Moreover, in the lattice momentum space the 

gauge condition (jD) becomes: 

p,Al = 0. (10) 
The most general structure of both diagonal and off-diagonal propagators is 

DM = (V - ^) D,{p') + ^D,{p') , (11) 

where Df i are the scalar functions. They are related to the formf actors D^jj{p) as follows: 

A(P^) = ^-f^D,^{p) , A(P^) = \{d,,{p) - D,{p')) . (12) 

It follows from (fTUI) that the longitudinal part of the propagator of the diagonal gluon, 
jjdi-aa ^ is zero. Thus we have three formf actors D^"^^ and D°^'^^^^. 




4 Numerical Results 

We calculate the propagators ((Tj), ((HI) on the symmetric lattices V = with L = 16, 24, 32 
using 50, 138 and 30 configurations, respectively. Simulations are done a.t (3 = 2.40 which 
corresponds to the lattice spacing a = (1.66 GeV)~^ j2I] if one fixes the physical scale 
^/a = 440 MeV. The details of the numerical gauge fixing procedure are given in the 
Appendix. 

We show all non-zero formfactors as functions of p"^ in Figure ^a) (in all Figures 
of this paper we depict the data obtained on 32^ lattice unless stated otherwise). One 
can see that all formfactors seem to tend to finite values as momentum goes to zero i.e. 
none of them is divergent or vanishing. The diagonal formfactor is dominating over the 
off-diagonal formfactors. 

In Figure Hfb) the gluon dressing function, p"^ D^p"^), is depicted. The diagonal (trans- 
verse) dressing function has a relatively narrow maximum at non-zero momentum Pg"^^ ~ 
0.7 GeV. Its behavior at small momenta is qualitatively very similar to the behavior of 
the gluon propagator in the Landau gauge (see e.g. [U|). The off-diagonal longitudinal 
dressing function has a wide maximum at p^^*^^^^ ~ 2 GeV, while for the transverse off- 
diagonal dressing function formfactor it is a monotonically rising function for all available 
momenta. 

One can compare the propagators obtained with the U{1) gauge conditions © and 
(jni). The first condition was implemented in Ref. fH] while the last one is adopted in the 
present paper. The comparison shows that the transverse diagonal propagators for these 
two gauge conditions coincide with each other at large momenta. At small momenta the 
formfactor obtained with Eq. (0) is slightly larger then the one calculated with Eq. ®. 
The difference at momentum p = 430 MeV is about 15%. The formfactors for the off- 
diagonal gluons coincide with each other for all available momenta. 

It is seen from Figure ^ that at j9 ~ 6 GeV dressing function D'^^^"-9,offdiag ^ g 
differ much from the free from 1/p^. The large renormalization of the formfactors in 
Landau gauge is discussed in Ref. [22] • There it has been shown that even three- loop 
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corrections do not describe the behavior of the lattice gluon propagator at p ~ 6 GeV. 

From Figures ma,b) it is clear that the off-diagonal gluon propagator is suppressed in 
comparison with the diagonal one. In Figure |21^a) we plot the ratio 

It is seen that the suppression of the off-diagonal propagator increases as the momentum 
decreases. This may be considered as an indication that in the MA gauge the diagonal 
gluons are responsible for physics in the infrared region. 

From Figuresn!^a,b) one may also notice that the formfactors Df^^^^^ip^) and D"^^^^^{p^) 
coincide at small momenta. In Figure Efb) we plot the ratio 



^offdiag(P^) - — offdiag. ox • (14) 



One can see that -Roffdiag decreases with decreasing momentum and vanishes at p ~ 1 
GeV. This implies that in the IR region the off-diagonal propagator has the form 

^^(p) ^ 5,. ■ Df'-^ip') , p'<lGeV. (15) 

In order to characterize the propagators quantitatively we have fitted the formfactors 
in the infrared region by the following functions: 

Dip')= ^itJj^^, (fit 2), (17) 

Dip')= (Yukawa fit), (18) 

Dip')= m-^^p4.pVm'^ > (Yukawa 2 fit) , (19) 

Dip')= (Gribovfit). (20) 

where Z, a , m and k are fitting parameters. The fitting functions (fT^ . (fT7|) . - after 
being modified to agree at large momenta with a known perturbation theory result - were 
used to fit the gluon propagator in Landau gauge [Hj. It was concluded that function (|T7jl 
provided the best fit. The fitting function (|T7jl was also used in the compact U(l) theory 
and compact Abelian Higgs model |2ni in three dimensions. 

The Yukawa fitting function, eq. (fTHj) . is introduced in order to compare our results 
with results of Refs. ^1 , where such behavior was assumed for off-diagonal propa- 
gators. Another interesting possibility is to consider the momentum dependent mass in 
eq. (|TH|l . — > m^(p^). Keeping only the lowest order of in m^(p^) expansion we get 
the fitting function (fTT^ where k is an additional dimensionless parameter. Finally, we fit 
the data by the function (plj) inspired by the Gribov proposal [211 the Landau gauge. 

The quality of the fitting result depends on the interval of momentum used in fitting. 
In this paper we restrict ourselves to the infrared region. We have chosen the interval 
starting from pmin = 27r/ (32a) = 0.325 GeV and ending at a variable value Pmax- For every 
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fit we determine Pmax as the highest momentum at which the data points corresponding 
to the lowest momenta are still consistent with the given fit. In other words, when 
momenta p > Pmax are included in the fit the fitting curves go off the error bars of the 
lowest momenta data points. We employ this procedure for the diagonal propagator only, 
because the statistical weight of the infrared data points is low^ while the significance of 
these points is high. In the case of the off-diagonal propagators we have limited the region 
of fit by highest momentum Pmax = 1-7 GeV. The fits are shown in Figures Efa-c) and 
the best fit parameters are presented in Tabled 

First we discuss the fits of the diagonal propagator. The three parameter fits (|16I17I19|) 
are working well and the corresponding curves are almost indistinguishable from each 
other. Figure Efa). However, the fitting function (fTTj) gives twice smaller value of Pmax 
than the other functions, indicating that fit ()17|) works in a narrower region than the other 
fits. The mass parameters for these fits do not coincide, and the difference between them 
is about 30%. 

We have also applied the two-parameter fits given by Yukawa ()18j) and Gribov pn|) 
functions. The Gribov fit is working well for the diagonal propagator. One can see from 
Figure El^a) that in order to discriminate the Gribov fitting function from the others we 
need the data at momenta smaller than available in our study. The Yukawa fit of the 
diagonal propagator does not work at all (we get x^/d-o.f. ~ 6 for fits in Pmax < 1 GeV 
region) . 

Concerning the transverse and longitudinal parts of the off-diagonal propagator one 
can make a few observations. First we notice, that the Gribov fitting function ()20|) is 
clearly not applicable for fitting of these propagators. Second, one can see that the 
formfactors for transverse and longitudinal parts almost coincide with each other at small 
momenta. The last fact implies that the best fit parameters for each particular type of 
the fits ()16II19|) must coincide as well, 

offdiag ofTdiag r7ofTdiaB ^ofFdiag ofTdiag offdiag offdiag offdiae 

^ ^ nil ' ~ ^/ 5 O^t ~ 5 l^t ^1^1 , 

^Therefore the x^-criterion can not be used for the definition of Pmax- 
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Figure 3: Fits of (a) transverse part of the diagonal propagator, (b) transverse and (c) 
longitudinal parts of the off-diagonal propagators, by the functions and (jl7|) . 



in agreement with Table (jT)). 

We have also found that the mass parameters in the off-diagonal gluon fits are ap- 
proximately two times bigger then the corresponding parameters in the diagonal fit prop- 
agators: 

Thus, the off-diagonal propagator is clearly short-ranged compared to the diagonal one. 

In Ref. ^3] the off-diagonal propagator was successfully fitted in the infrared region 
of coordinate space by the Yukawa propagator. Our results (with respect to the off- 
diagonal propagator) indicate, that other fitting functions can also be used to describe 
this propagator. 

Summarizing, we can say that the diagonal propagator can be fitted almost equally 
well with any of the three-parameter fits ()16p. ()17|) and ()19|). The same is true for the 
off-diagonal propagator. Among the two-parameter fits ()20j) is superior for the diagonal 
propagator while ()18|) is better for the off-diagonal one. 

At smallest available momentum (325 MeV) the off-diagonal propagator is suppressed 
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fit 


m, GeV 


a OT K, 


Z 


Pmax, GeV 


x'/dof 


Transverse diagonal 


fit 1 


0.73(2) 


0.92(3) 


16.9(4) 


1.7 


0.8 


fit 2 


0.58(2) 


0.49(5) 


8.5(2) 


1.0 


0.4 


Gribov fit 


0.33(1) 




4.58(5) 


0.9 


0.9 


Yukawa 2 fit 


0.50(2) 


0.19(3) 


8.3(3) 


1.7 


0.9 


Transverse off- diagonal 


fit 1 


1.6(2) 


0.6(2) 


1.3(2) 


1.7 


1.0 


fit 2 


1.26(4) 


0.19(4) 


0.73(2) 


1.7 


1.0 


Yukawa fit 


1.08(2) 





0.63(1) 


1.7 


1.5 


Yukawa 2 fit 


1.29(6) 


0.15(5) 


0.81(5) 


1.7 


1.0 


Longitudinal off-diagonal 


fit 1 


1.4(2) 


0.5(3) 


1.0(3) 


1.7 


1.1 


fit 2 


1.14(6) 


0.19(6) 


0.63(3) 


1.7 


1.1 


Yukawa fit 


0.96(3) 





0.52(1) 


1.7 


1.3 


Yukawa 2 fit 


1.14(8) 


0.12(7) 


0.66(6) 


1.7 


1.1 



Table 1: The fitting results for the propagator formfactors at low momenta: best param- 
eters of the fits (fTT)H2IHl . The corresponding highest momentum, Pmax, and x^/d-o.f. are 
also presented. 

by the factor about 50 with respect to the diagonal one. Note that this suppression 
is only partially due to larger values of the mass parameter m while equally or even 
more important role is played by smaller values of the parameter Z. We expect that 
the similar suppression exists also in the continuum theory because at our lattice spacing 
the renormalization effects for the Z-parameter are already quite small (according to 
Figure H^a) the transverse diagonal and off-diagonal formfactors almost coincide with 
each other at largest available momentum). 

5 Conclusions 

Our results obtained in the Maximally Abelian gauge of SU{2) gluodynamics clearly 
show that at low momenta the propagator of the diagonal gluon is much larger than the 
propagator of the off-diagonal gluons. This suggests that the colored objects at large 
distances interact mainly due to exchange by the diagonal gluons in agreement with the 
Abelian dominance property established in numerical studies of the MA gauge ^Hj for 
fundamental test charges^. 

The propagators do not show indications that they are either vanishing or divergent 
when p ^ 0. To provide a quantitative description of the propagators at low momenta we 
fit the propagator formfactors using various functions. All infrared fits for both diagonal 
and off-diagonal propagators contain massive parameters, which are non-zero for both 

•^As for adjoint charges see discussion in Ref. |25| . 
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propagators. When the same fitting function is applied to the diagonal and off-diagonal 
formfactors the mass parameter for off-diagonal formfactor is more than twice bigger than 
that for the diagonal one. This is in a qualitative agreement with findings of Ref. [14 . But 
our more detailed analysis revealed that the difference in values of the mass parameters is 
not the only reason of the off-diagonal propagator suppression, the other reason is small 
values of parameter Z . 

We have found that the diagonal propagator has qualitatively the same momentum 
dependence as the gluon propagator in Landau gauge while the off-diagonal propagator 
is very different. At small momenta the off-diagonal propagator is diagonal with respect 
to the space indices and thus defined by a single scalar function since its transverse and 
longitudinal formfactors become equal within error bars. 
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Appendix 

The Maximally Abelian and the Abelian Landau gauges are defined as global maxima of 
the functionals Q and respectively. The global maxima is difficult to reach numeri- 
cally and usually one finds several field configurations corresponding to local maxima of 
the gauge fixing functional and then the configuration with the highest value of the func- 
tional is chosen. The choice of the correct maximum, known as a Gribov problem (25, is 
crucial for the observables both in the MA gauge of SU{2) gauge model j2E] and in the 
Landau gauge of the U(l) gauge model P7j . 

Gauge fixing of the MA gauge with a careful treatment of the Gribov ambiguity 
was studied in Ref. [221 • In our investigation we follow this paper using the Simulated 
Annealing algorithm with 10 randomly generated gauge copies. The MA gauge fixing 
procedure is described in j2ni, and the Abelian Landau gauge fixing algorithm is briefiy 
considered below. 

To maximize the functional © we use the local over-relaxation algorithm with u = 1.8, 
see, e.g., Ref. with 20 randomly generated gauge copies. For the local gauge fixing 
we choose the following convergence criterion: 

4 

Gi^) = I cosLPf^{x) sm9f^{x) - cosLp^{x + fx) sm9^{x + < e , (Al) 

where e is a small parameter. Note, that this condition must be satisfied at each site x 
of the lattice. 
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Figure 4: (a) The longitudinal part of the propagator of the diagonal gluon at various 
momenta vs. convergence parameter (b) the volume dependence of the diagonal and 
off-diagonal propagators at p = 0, the corresponding lattice volumes are: 16^, 24^, 28^ 
and 321 



To study the effects of the incomplete gauge fixing we chose the longitudinal part of 
the propagator of the diagonal gluon D'l^^^ because it is most sensitive to the details of the 
gauge fixing procedure. According to the local gauge condition (jH), D^^^^ must be zero 
when perfect numerical procedure is used. Its dependence on the convergence parameter 
e is presented in Figure |3fa). The propagator significantly depends on e, especially in the 
region of small momenta. In our simulations we choose e = 10~®. 

We have also checked the dependence of Df^'^^ on the number of gauge copies, (Ngc), 
used both in the MA and in the Abelian Landau gauge fixings. This check has been done 
using set of 30 configurations on L = 24 lattice. Within error bars we have observed no 
dependence on the number of the MA gauge Gribov copies and we have found only very 
mild dependence on the number of the Abelian Landau gauge copies. 

To check the finite volume effects we have calculated the propagators on different lat- 
tices at the same /3 = 2.40. The transverse formfactors at zero momentum /}^*'^9:o//'^«agj^g-j 
were calculated as follows: 

£fdia.g,offdiag _ }^ jjfdiag,of fdiag {A.2) 

We find that within error bars the values of the zero-momentum propagator are indepen- 
dent of the lattice volume, as can be seen from Figure lU^b). 

We have also studied the volume dependence of the transverse part of the diagonal 
propagator at non-zero momenta. Since the finite volume affects mainly the low momen- 
tum region we have concentrated on the propagator at p < 1.5 GeV. In Figure El we plot 
the transverse part of the diagonal propagator calculated on 24'^ and 32'^ lattices. It is 
seen that the volume dependence is very weak, it is within the error bars. The position 
of the maximum of p^D{p'^) is the same for both lattices. The fit of the data for the 
propagator on 24^ lattice gives parameters which differs less than by 4% from that given 
in Table 1 for 32^ lattice. Thus we estimate the systematic error induced by the finite 
volume effects to be less than 4%. 
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Figure 5: The transverse diagonal propagator at small momenta for the lattices L = 24 
and L = 32. 
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